Slow dynamics in a 2D Ising model with competing interactions 
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The far-from-equilibrium low-temperature dynamics of ultra-thin magnetic films is analyzed by 
using Monte Carlo numerical simulations on a two dimensional Ising model with competing exchange 
(Jo) and dipolar (Jd) interactions. In particular, we focus our attention on the low temperature 
region of the (S, T) phase diagram (where S — Jo/Jd) for the range of values of 6 where striped 
phases with widths h = 1 (hi) and h — 2 (h2) are present. The presence of metastable states of 
the phase h2 in the region where the phase hi is the thermodynamically stable one and viceversa 
was established recently. In this work we show that the presence of these metastable states appears 
as a blocking mechanism that slows the dynamics of magnetic domains growth when the system is 
quenched from a high temperature state to a low temperature state in the region of metastability. 

PACS numbers: PACS numbers: 75.40.Gb, 75.40.Mg, 75.10.Hk 
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I. INTRODUCTION 

In the last years a strong effort has been devoted to 
understand both the equilibrium and out-of equilibrium 
properties of ultrathin magnetic filmsi. These materials 
have attracted much attention mainly due to their po- 
tential applications, such as information storage 2 . Ultra- 
thin films find also many important applications both in 
biotechnology and pharmacology. It is today a well estab- 
lished experimental fact that the magnetization processes 
in ultra-thin magnetic films are ruled by the microscopic 
competition between short range ferromagnetic couplings 
and long-range frustrated antiferromagnetic dipolar in- 
teractions, which give place to very novel dynamical and 
static behaviorsi. 

It is worth mentioning that both the theoretical and 
the experimental interest in studying systems with com- 
petition between short-range ordering interactions and 
long-range frustrating interactions widely exceeds the 
field of ultra-thin films. Actually, many different ex- 
perimental systems can be modeled by this kind of mi- 
croscopic interactions, which give place to very rich dy- 
namical and static properties. In soft-matter physics for 
instance, we can mention diblock copolymer melt and 
cross-linked polymer mixtures, among others. Type I 
superconductors and rare-earth layers that occur in the 
perovskite structure of REBa2Cu3C>7 — 6 (where RE rep- 
resents a rare earth from the lanthanidc series) can be 
very well modeled with these interactions^. 

It has also been frequently suggested that competing 
interactions can explain many of the phenomenological 
features observed in the glass formation process and in 
supercooled liquids^. Summarizing, many of the conclu- 
sions drawn from this work can be surely be applied to a 



large variety of physical systems. 

For sufficiently thin films the magnetic moments align 
perpendicular to the plane of the film, indicating that the 
surface anisotropy is sufficient to overcome the anisotropy 
of the dipolar interaction which favors in-plane ordering. 

Works in two dimensional uniaxial spin systems, where 
the spins are oriented perpendicular to the lattice and 
coupled with these kind of interactions, have shown a 
very rich phenomenological scenario concerning both its 
equilibrium statistical mechanics^ and non-equilibrium 
dynamical properties^i In particular, some of these 
results 2 '- showed the existence of different types of slow 
relaxation dynamics when the system is quenched from a 
disordered high temperature configuration to a subcrit- 
ical temperature, depending on the relative strengths of 
the dipolar and exchange interactions. 

Under these circumstance one can use a uniaxial Ising 
representation for describing the magnetic moments^. 
The ultra-thin film is then described by the Hamiltonian 
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where the spin variable Cj = ±1 is located at site i of 
a square lattice, the sum Yluj) runs over all pairs of 
nearest neighbor sites and Y)u j-\ runs over all distinct 
pair of sites of the lattice; is the distance (in crystal 
units) between sites i and j, 8 represents the quotient 
between the exchange Jo and dipolar Jd coupling pa- 
rameters (8 = Jo/Jd). The energy is measured in units 
of Jd, which is always assumed to be antiferromagnetic 
(Jd > 0). Hence 8 > means ferromagnetic exchange 
coupling. 



We have recently studied in detaiUS the low tempera- 
ture phase diagram of this system in the region where the 
change in the relaxation properties has been observed. 
We showed that for very low temperatures metastable 
states appear. In this work we investigate the effects 
of the presence of these metastable states on the far- 
from equilibrium dynamical properties of the system. In 
section [H] we present a review of the equilibrium phase 
diagram and metastability properties in the region of in- 
terest. In section IIIII we analyze the magnetic domain 
growth or coarsening dynamics of the system when it 
is quenched from a disordered state (which corresponds 
to infinite temperature) to a temperature below the or- 
dering transition for different values of 8. Using Monte 
Carlo simulations we study the statistics of domains of 
the striped phases hi and h2. We then analyze the tem- 
poral behavior of the average linear size of the domains 
L. We show that the coarsening dynamics is strongly 
affected by the presence of metastable states, which gen- 
erate blocking clusters of the metastable phase where the 
domain walls of the stable phase become pinned. Such 
blocking clusters generate free-energy barriers to the do- 
main growth dynamics that are independent of the linear 
domain size. Some conclusions and remarks are summa- 
rized in section Hvl 
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FIG. 1: Phase diagram (5, T) in the region of parameters 
under study. Filled triangles correspond to the critical tem- 
peratures T c (8) obtained by specific heat calculations for the 
phase transition between the ordered antiferromagnetic (AF) 
and striped phases hi and h2 and the paramagnetic (para) 
one. Filled circles (open squares) correspond to the stability 
line of the hi (h2) phase, obtained by analyzing the staggered 
magnetization Mhi (Mh2)- Filled diamonds correspond to the 
first order transition lines between the hi and h2 phases and 
also between the hi and AF phases, obtained by the free en- 
ergy numerical calculations. The shaded region indicates the 
presence of metastable states. 



II. EQUILIBRIUM PHASE DIAGRAM AND 
METASTABLE STATES 

The overall features of the finite temperature phase di- 
agram associated with Hamiltonian Q were described by 
Maclsaac and coauthors^ by means of Monte Carlo sim- 
ulations on 16 x 16 lattices and analytic calculations of 
the ground statei. They found that the ground state 
of Hamiltonian Q is the antiferromagnetic state for 
8 < 0.425ii. For <5 > 0.425 the antiferromagnetic state 
becomes unstable with respect to the formation of striped 
domains structures, that is, to state configurations with 
spins aligned along a particular axis forming ferromag- 
netic stripes of constant width h, so that spins in adja- 
cent stripes are anti-aligned, forming a super lattice in 
the direction perpendicular to the stripes. At high tem- 
peratures, of course, the system always becomes param- 
agnetic. Specific heat calculations showed that the tran- 
sition between the paramagnetic and the striped phases 
is a second order onei 

We have recentlj*^ performed Monte Carlo simulations 
of Hamiltonian on square lattices up to 48 x 48 sites 
using periodic boundary conditions and heat bath dy- 
namics. Our calculations focused on the low tempera- 
ture region of the (8, T) phase diagram for values of 8 
between 0.2 and 2, which includes the transition line be- 
tween the AF and the striped phase with width h = 1 
(hi) and also the transition line between the hi phase 
and the striped phase with width h = 2 (h2). First we 
calculated, through the energy fluctuations, the specific 
heat C as a function of temperature for different values 



of 8 and different system sizes up to 48 x 48 sites. By 
considering the peaks in the specific heat we obtained the 
second order critical line between the paramagnetic and 
the low temperature ordered phases hi and h2. These 
results (see Fig. slightly improved those obtained by 
Maclsaac and coauthors 3 for 16 x 16 lattices, thus show- 
ing a fast convergence of the critical temperature for in- 
creasing system sizes, at least for small values of 8. 

Next, through a numerical study of the free energy, we 
analyzed the transition between the hi and h2 phases. 
The free energy of the phase hi (h2) was calculated for 
increasing (decreasing) values of 8. We observed a contin- 
uous change of the minimal free energy from one phase 
to the other, with a discontinuous change in the slope 
for 8 — 1.26(1) indicating the presence of a first order 
phase transition 10 . This transition is indicated by means 
of diamonds in Fig. ^ We have also repeated these cal- 
culations for the transition line between the hi and AF 
phase, finding similar results. 

Close to 8 — 1.26(1) the free energy displays a multi- 
valued behavior characteristic of a first order phase tran- 
sition. This behavior signals the metastable nature of 
these phases in some parts of the phase diagram. To 
characterize the presence of metastable states observed 
in the transition between the hi and h2 phases we intro- 
duced the staggered magnetizations Mhi and M^2, and 
also their associated susceptibilities Xhi and Xhi for the 
hi and h2 phases^ -. These quantities permitted us to 
analyze the stability of both phases in the different parts 
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of the phase diagram. In Fig. Q]the shaded region indi- 
cates the presence of metastable states. It is important to 
stress that in the shaded region inside phase hi the only 
phase observed to be metastable was hi. For 6 > 1.26 
metastables states of phases of higher width were ob- 
served, leading to a much more complicated metastable 
region. In this work we will focus only on the dynami- 
cal behavior of the system in the hi region. As we will 
show in the next section, for a fixed value of 6 inside 
this region different dynamical regimes are observed as 
the temperature is lowered and one enters the region of 
metastability. 
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III. SLOW DYNAMICS 

When a system is quenched from a high temperature 
disordered phase into a low temperature ordered phase 
domains form and grow, a process that is known as coars- 
ening. The coarsening process has been extensively stud- 
ied both experimentally and theoretically over the past 
decade^. 

Perhaps the most thoroughly studied system, and also 
the most common example is the Ising model. When this 
system is quenched from a high temperature to one below 
its critical temperature (T < T c ) ferromagnetic domains 
of up and down spins form and coarsen. The system 
presents curvature driven growth and the characteristic 
domain size L grows with time as L(t) ~ t 1 / 2 . If the 
system is cooled to zero temperature the domain walls 
can be easily determined as bonds between oppositely 
oriented spins, but if the system is cooled to a temper- 
ature different from zero it becomes difficult to define 
domains and domain walls since small islands generated 
by thermal fluctuations arise. To overcome this problem, 
Derridai£ proposed a new method to measure properties 
related to coarsening in the presence of thermal fluctu- 
ations. This method was extended by Hinrichscn and 
Anton:— to determine domain walls for nonzero temper- 
atures. The method compares the state of a system with 
replicas in the different ground state configurations when 
they are all submitted to the same thermal noise, that 
is, when the same sequence of random numbers is used 
to update all systems. In this way, if one starts from 
a replica in the ordered state a spin flip will be a con- 
sequence of the thermal noise. When a spin flip occurs 
simultaneously in all the replicas it can be considered as 
a thermal fluctuation, otherwise the fluctuation will be 
due to the coarsening process. We used this technique 
to study the dynamics of domain walls, and character- 
ize the coarsening process, when the system described 
by Hamiltonian Q is quenched from a high temperature 
disordered state into the region where it orders. In par- 
ticular we focus our interest on the growth of domains 
of the striped phase hi when the metastability line is 
crossed for values of 0.8 < 8 < 1.26 (See FigQJ. 

To characterize the growth of the domains we de- 
termined first the domain areas A(t), by counting the 



FIG. 2: Characteristic domain size L(t) vs. t when 5 — 1.1 
and T = 0.2 for three different system sizes N = 24 x 24, 
N = 36 x 36 and N = 48 x 48 . The dashed line indicates 
L(t) ~ i 1//2 . The inset shows the data collapse obtained using 
finite size scaling analysis. 



number of spins inside each domain. The characteristic 
(linear) domain size was calculated as L(t) = \J (A(t)), 
where (A(t)} is the mean domain area of the system at 
time t. 

In Fig. |21 we present the behavior of the characteristic 
domain size L(t) when 8 = 1.1 and T = 0.2 for three 
different system sizes N = 24 x 24, N = 36 x 36 and 
N = 48 x 48. After a short transient in which the char- 
acteristic length presents a slow growth the system enters 
into a coarsening regime where L(t) ~ i 1 / 2 , as expected 
for a system with non-conserved order parameter—. For 
large times L(t) presents a crossover to a saturation 
value. This saturation behavior is clearly a finite size 
effect, since a domain cannot grow beyond the system 
size N so that as the system size increase, the behavior 
L(t) ~ i 1 / 2 remains for larger periods of time. 

In section [H] we observed and characterized the pres- 
ence of metastable states in the low temperature region 
of the phase diagram. We will study now how the behav- 
ior of L(t) changes as we lower the temperature and cross 
the metastability line for a fixed value of 8. In Fig. [3] we 
present the time evolution of L(t) when 8 = 1.1 for eight 
decreasing temperatures. For low temperatures a regime 
of slow growth develops at intermediate time scales be- 
fore the system crosses over to the t 1 / 2 coarsening regime. 
As we lower the temperature the intermediate regime ex- 
tends to larger time scales, but all the curves eventually 
cross over to the L(t) ~ i 1 / 2 regime. Note that for short 
times there seems to be a change in the concavity of L(t) 
when it crosses T — 0.1, which coincides with the bound- 
ary of the metastable phase for 8 = 1.1. 

To characterize the behavior of L(t) as we lower the 
temperature we studied the crossover time r from the 
power law to the saturation regime. Note that the 
crossover time to the power law regime presents a similar 
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FIG. 3: Characteristic domain size L(t) vs. t when 5 = 1.1 
(N — 24 x 24) for eight different temperatures, starting form 
the left curve T = 2.0, 0.175, 0.15, 0.1, 0.085, 0.075, 0.065 
and T = 0.05. Note that all the curves eventually cross-over 
to the L(t) ~ t 1 / 2 regime indicated with a dashed line. 



behavior. However the study of the intersection of the 
power law branch of the curve and the horizontal satu- 
ration branch allows for a sistematic approach. This is 
so since the saturation value always correspond to the 
linear size of the system. In figure 21 we present how r 
grows as the temperature is lowered. For temperatures 
greater than T = 0.1 the crossover time presents a linear 
dependency with 1/T, while for temperatures lower than 
T = 0.1 it presents an exponential increase with 1/T as 
can be observed in the Arrhenius plot presented in the 
inset of Fig. ^ 

The straight line indicates the best fit, given by a func- 
tion of the form 



t = T exp(ri/T) 



(2) 



where To = 62.5(5) and t\ = 0.39(5). Using this expres- 
sion we present in Fig. [3 a data collapse plot of Lit) in 
the low temperature regime. 

In the high temperature regime the crossover time de- 
creases linearly as the temperature increases, and L(t) 
collapses simply by scaling with T as can be seen in fig- 
ure 

Summarizing, for 5 = 1.1 two different dynamical 
regimes were observed above and below the metastability 
line (T = 0.1). When the system is quenched to the or- 
dered phase to a temperature T > 0.1, the characteristic 
domain size L(t) grows as i 1 / 2 after a short transient. If, 
on the other hand, the system is quenched to a tempera- 
ture T < 0.1, an intermediate regime with a slow growth 
appears. For long times the system always crosses over 
to the t 1 / 2 regime. That is, for every temperature we ob- 
serve the same behavior presented in FigO where L(t) al- 
ways reaches the asymptotic behavior t 1 ! 2 as we increase 
the system size. However, as we lower the temperature 
the crossover time to this regime increases. We repeated 
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FIG. 4: Crossover time r vs T for 8 = 1.1 and N = 24 x 24. 
On the inset an Arrhenius plot r vs 1/T for the four lower 
temperatures is presented. The straight line indicates the best 
fit. 



these analysis for different values of 5, inside the h = 1 
region, obtaining the same qualitative behaviors. 

These dynamical behaviors present a strong resem- 
blance with the ones observed in the two dimensional 
Shore modeli^. This model is a ferromagnetic Ising 
model on a square lattice with frustration added by in- 
troducing weak next-nearest-neighbor antiferromagnetic 
bonds, that is 
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(3) 



The presence of NNN antiferromagnetic bonds in this 
model introduce free-energy barriers to domain coars- 
ening that are independent of the domain size Li^. 
Such barriers in this model are a consequence of a cor- 
ner rounding process which generates structures that 
block the coarsening dynamics^. Hence, the system is 
stuck and coarsens little on time scales t <C tb(T) — 
exp(FB /T) (Fb being the height of the barrier), while 
on time scales t 3> tb (T) the free-energy barrier can be 
crossed and the t 1 / 2 behavior emerges. 

Thus, it is interesting to see if structures that block the 
coarsening process are present in our system. In fact the 
presence of metastable states for low temperatures gives 
a strong hint of what to look for. 

Figure presents a series of snapshots of the coarsen- 
ing process when a system of size N = 48 x 48 and 5 = 1.1 
is quenched to T = 0.05. For these values of 5 and T the 
system evolves in the metastability region. Two different 
kind of domain walls can be clearly distinguished. When 
domains of phase hi with different orientation meet, the 
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FIG. 5: Data collapse plot of L(t) in the low temperature 
regime for three different temperatures (8 = 1.1 and N = 
24 x 24). The crossover time r corresponds to the best fit 
presented in Eq. JSJ. 



spins in the h2 phase. This coarsening behavior presents 
a strong resemblance to the one observed in the snapshots 
presented in Fig. [7| 

Figure El shows the time r/ l2 to shrink a 4 x 4 square 
in phase h2 immersed in a system with N = 48 x 48 
spins in phase hi with fixed boundary conditions. As 
the temperature is lowered the shrinking time diverges 
as Th2 ~ exp(0Al/T). This divergence agrees well with 
the divergence observed in the crossover time from the 
slow growth to the i 1 / 2 regime. We have also studied the 
shrinking time of blocks in different phases, such as com- 
binations of vertical and horizontal hi phases and also 
a ferromagnetic block immersed in the hi phase. In all 
these cases a similar behavior was observed. The sys- 
tem quickly reached a configuration where small blocks 
of the h2 phase were present, with a slowing down of the 
domain growth similar to the one described above. 



IV. CONCLUSIONS 
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FIG. 6: Data collapse plot of L(t) in the high temperature 
regime for three different temperatures S — 1.1 and N = 
24 x 24. 



walls detected by the method are thin. The thicker walls 
correspond to bands of two spins, that is stripes of phase 
h2. We should stress that phase h2 is the only metastable 
phase in phase hi—. Notice how these small blocks of the 
h2 phase seem to slow the domain growth by pinning the 
domain walls. These small blocks were also observed at 
short times when the system was quenched to a tempera- 
ture above the metastability line, however they are highly 
unstable and did not seem to block the domain growth. 
In order to quantify this effect, we studied the time for 
completely shrinking squares in phase h2 immersed in 
a system in phase hi with fixed boundary conditions. 
A square was considered to be completely shrunk when 
the length of the sorrounding domain walls became zero. 
FigJSJ presents a series of snapshots of this process. In 
the first time steps the square quickly deforms. However, 
further advance seems to be blocked by small blocks of 



In this paper we have studied some dynamical prop- 
erties of a two-dimensional Ising Hamiltonian with com- 
peting interactions. We analyzed the coarsening process 
when the system is quenched from a high temperature 
disordered phase into the ordered phase hi. To char- 
acterize the growth of domains we considered the time 
evolution of the characteristic domain size L(t). We 
found that for a fixed value of S the system presents 
two different dynamical behaviors associated with the 
presence or absence of h2 metastable states. When the 
system is quenched into the ordered phase for temper- 
atures above the metastability region, the characteristic 
domain length presents a power law t 1 / 2 growth. If, on 
the other hand, the system is quenched to temperature 
in the metastability region, the behavior of Lit) presents 
a slow growth intermediate regime before crossing over 
to the i 1 / 2 power law growth. Through a direct exam- 
ination of snapshots of the system during the coarsen- 
ing process we found that the presence of small domains 
in the h2 phase slowed the coarsening when the system 
was quenched to the metastable region. These results 
are consistent with the presence of free-energy barriers 
independent of the domain size L, associated with block- 
ing clusters of the metastable phase, which generates a 
crossover in the coarsening behavior as we cross the spin- 
odal line. Since the cross over time diverges as the tem- 
perature is lowered, the very slow behavior at interme- 
diate times may be indistinguishable from a logarithmic 
law. This could explain the apparently logarithmic scal- 
ing observed in the aging behavior of this model^. 
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sejo Nacional de Investigaciones Cienti'ficas y Tecnicas 
CONICET (Argentina), Consejo Provincial de Investiga- 
ciones Cienti'ficas y Tecnologicas (Cordoba, Argentina), 
Secretaria de Ciencia y Tecnologia de la Universidad Na- 
cional de Cordoba (Argentina) and Fundacion Antorchas 
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FIG. 7: Snapshots of the coarsening process for S = 1.1, N — 48 x 48 and T = 0.05. The squares at the top correspond to 
spin configurations at different times of the coarsening (black points: up spins; white points: down spins). The black areas in 
the figures at the bottom are the corresponding domains of the hi phase at every time, while the white lines correspond to the 
domain walls. 



(Argentina). 
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FIG. 8: Snapshots of the coarsening process for a 4 x 4 shrink- 
ing square in phase hi immersed in a system in phase hi with 
fixed boundary conditions. 
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FIG. 9: Shrinking time of a 4 x 4 square in phase hi immersed 
in a system in phase hi. 
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